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ABSTRACT 

We have used the Australia Telescope Compact Array (ATCA) and the Swedish-ESO Sub- 
millimetre Telescope (SEST) to map the large-scale atomic and molecular gas in the nearby 
(4 Mpc) Circinus galaxy. The ATCA H I mosaic of Circinus e xhibits the warps i n position 
angle and inclination revealed in the single-pointing image of iJones et alJ dl999l) . both of 
which appear to settle beyond the inner 30 kpc which was previously imaged. The molecular 
gas has been mapped in both the CO J = 1 — » and J = 2 — > 1 transitions down to 
a column density of Nn 2 > 10 21 cm" 2 (3er), where we derive a total molecular gas mass 
of Mh 2 ~ 2 x 10 9 M Q . Within a radius of 3 kpc, i.e. where CO was clearly detected, the 
molecular fraction climbs steeply from w 0.7 to unity (where Nn 2 — 4 x 10 22 cm" 2 , cf. 
iVm = 10 21 cm -2 ) with proximity to the nucleus. Our Hi mosaic gives an atomic gas mass 
of A/hi ~ 6 x 10 9 M Q , which is 70% of the fully mapped single dish value. Combining 
the atomic and molecular gas masses gives a total gas mass of M gas = Mm + Mh 2 ~ 
1 x 10 10 M Q , cf. the total dynamical mass of rs 3 x 10 11 M Q within the inner 50 kpc of 
our mosaiced image. The total neutral gas mass to dynamical mass ratio is therefore 3%, 
consistent with the SAS3 classification of Circinus. The hi gh (molecular) gas mass fraction 
of Mh 2 /Mdyn ~ 50% found previously dCurran et alJI 19981) . only occurs close to the central 
w 0.5 kpc and falls to < 10% within and outwith this region, allaying previous concerns 
regarding the validity of applying the Galactic -/Vh 2 /7co conversion ratio to Circinus. The 
rotation curve, as traced by both the Hi and CO, exhibits a steep dip at « 1 kpc, the edge of 
the atomic/molecular ring, within which the star-burst is occurring. We find the atomic and 
molecular gases to trace different kinematical features and believe that the fastest part (> 130 
kms -1 ) of the sub-kpcring consists overwhelmingly of molecular gas. Beyond the inner kpc, 
the velocity climbs to settle into a solid body rotation of « 150 kms" 1 at > 10 kpc. Most 
of the starlight emanates from within this radius and so much of the dynamical mass, which 
remains climbing to the limit of our data (> 50 kpc), must be due to the dark matter halo. 

Key words: galaxies: active — galaxies: individual: Circinus — galaxies: kinematics and 
dynamics — galaxies: ISM — galaxies: Seyfert — radio lines: galaxies. 



1 INTRODUCTION 

Circinus is a nearby, highly inclined spiral galaxy located behind 
the Galactic Plane (at I, b = 311?3, -3? 8), whe re the high stellar 
density and dust extinction (j4b = 6.3 mag, ISchlegel et alJI 19981) 
make it difficult to study its optical properties and enviro nment. At 
a distance of only 4.2 ± 0.8 Mpc JFreeman et alJI 19771) . Circinus 
represents one of the closest examples of a star-burst galaxy with 
a Seyfert2 nucleus (Moorwood & Oliva 1990; Oliva et al. 1994). 
Low-resolution H I observations of Circinus with the 64-m Parkes 
radio telescope revealed a vast galaxy with H I 21 -cm emission ex- 
tending to a half-width of 32' JFreeman et alJI 19771) . Its enormous 
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extent has also been noted in both the shallow and deep H I sur- 
veys of the Zone of Avoida nce (HIZSS, Henning et al. 2000 and 
HIZOA; ljuraszek et al J 1 20001 respectively) and, most recently , in 
the Hi Parkes All-Sky Survey (HIPASS. [Koribalski et alj|2004l) . 

Hig h resolution, sing le-pointing ATCA Hi observations of 
Circinus Uones et al J 19991) reveal a complex velocity field and gas 
distribution out to a diameter of ~ 80' (~ 100 kpc), with strong 
kinematical warping (also evident between the large-scale position 
angles and inclinations, Table[T](. The atomic gas is distributed in a 
~10 kpc ra dius ring enclosi ng a bar, in turn terminating in a 1 kpc 
radius ring llones et alJI 19991) . This inner H I ring may be an outer 
component of the 600 pc C O ring dCurran e t al. 1998), where the 
star formation is dominant (Mar coni et al.lll994T) . This itself may 
feed the nucleus via a nuclear bar nestled within the molecular ring 



2 S. J. Curran, B. S. Koribalski & I. Bains 



Table 1, General properties of the Circinus galaxy. 



central position 


"J2000 


14 h 13 m 09!95 (±0?02) 


1 




"5.J2000 


-65°20'21'.'2(±0'.'l) 




type 




SAS3 


2 


adopted distance 


D 


4.2 Mpc 


3 


optical extent* 


D24.5 


11.'9 x 4f6 


3 




£>26.6 


17.'2±lf7 


3 


position angle 


PA opt 


220° 


2 




PAco 


214° ± 4° 


4 




PA m 


210° ± 5° 


3 


inclination 


^opt 


64° 


2 




*co 


78° ± 1° « 30") 


4 




'III 


65° ± 2° 


3 


systemic velocity 


''III 


439 ±2 km s" 1 


3 


H I flux density 


F m 


1910±130Jykms- 1 


3 


H I mass 


M m 


7.9 ±0.5 x 10 9 M Q 


3 


dynamical mass 


Afdyn 


1.3 ±0.2 x 10 11 M© 


3 


B-band luminosity 


L B 


7 ± 2 x 10 9 L 


3 


B-band extinction 


A B 


6.3 mag 


5 



Notes: Freeman et al. (1977) derived extinction-corrected optical diameters 
for the Circinus galaxy at the 24.5th blue magnitude, D24.5, and, through 
extrapolation, the 26.6th blue magnitude, D26.6 (the Holmberg diameter). 
References: ( 1 ) Greenhill et al. 1 2003^ (2)lde Vaucouleurs et aljjl99ll), (3) 
1 et all il9' 



iFreeman et al] Jl977l) . (4) ICurran et a l. (1998), (5) Schlegel et al. 1 1998). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 H 1 observations 

The H I observations were performed with the ATCjAy in its 375-m 
configuration, on June 18-22, 1999. In order to avoid solar inter- 
ference on the short baselines of the array, which are needed to 
detect the very extended H I emission, we only observed at night. 
We mosaiced Circinus using a grid of 5 x 3 points along the major 
axis with pointing centres separated by 16'5 (i.e. half the primary 
beam-width) in order to fully Nyquist sample the field. The band 
was centered at 1418 MHz with a bandwidth of 8 MHz, divided 
over 512 channels, giving a channel spacing of 3.3 kms -1 , i.e. a 
velocity resolution of 4 kms -1 . We observed for 5 x 12 hours in 
total, achieving an integration time of 44 hours on Circinus, i.e. 
~3 hours per pointing. We used PKS B 1934— 63 for flux and band- 
pass calibration and PKS B 1329-665 for phase calibration. Data 
reduction and analysis were performed with the MIRIAD software 
package. 

After identifying the velocity range of H I emission in Circi- 
nus and noting the foreground emission from the Galaxy and high 
velocity clouds, we subtracted an average of the line-free channels 
from the wv-data. The H I data were then Fourier-transformed us- 
ing 'natural' weighting. We used ten baselines between 31m and 
459 m, excluding the five longer baselines to the distant antenna six, 
which do not contribute to the H I emission. The task mossdi was 
used to clean the mosaiced H I cube which was then restored with 
a beam size of 124" x 107". The r.m.s. noise per each 4 kms~ 
channel is ~ 2 mjy beam -1 , close to the theoretical value. 



JMaiolino et al.ll2000f) . thus providing a means of transporting gas 
from the outer galaxy to the central engine (e.g. IShlosman et al] 
119891) . 

The notion of a continuous structure from the kpc to sub-pc 
scale is s upported by the VLB I observations of the H2O maser 
emission (Greenhill et al. 2003) and VLTI near infrared observa- 
tions ((Tristram et al. 2007), which both show that the sub-pc accre- 
tion disk/obscuring torus shares a simi lar close to edge-on orienta- 
tion to the larger-scale molecular gas JCurran et alj|l998l) . The in- 
ner edge of the warped maser disk shares a similar position angle to 
the large-scale gas disk (PA = 209° ±3 °), although this increases 
to 236° ± 6° at the outer edge (Greenhill et al. 2003). Furthermore, 
a second population of masers are found to trace a wide angle out- 
flow which coincides with the > 100 pc-scale ionisa tion cone (Mar- 
coni et a l. 1994, Veilleux & Bland-Hawthornlll997t) and molecular 
outflow (Curran et al. 1999), as well as sharing the same position 
angle as the radio lobes (Ha rnett et alJl99ulElmouttie et al.l 19951) . 
The lobes, cone and outflow are all directed along the minor axis of 
the galaxy, coincident with the ro tation axis of the molecular ring 
dCurran et a0l999l;ICurranll200lh. 

Previously, I Jones et al] J19991) obtained very sensitive, single- 
pointing H I data of the Circinus galaxy with the Australia Tele- 
scope Compact Array (ATCA). Since the ATCA primary beam is 
~33', these data were not sufficiently sensitive to detect the outer 
edges of the large gaseous disk of Circinus. Here we present and 
analyse an H I mosaic of Circinus, obtained with the 375-m config- 
uration of the ATCA, covering the whole galaxy. This is comple- 
mented by the most extensive CO 1 — ► and 2 — » 1 maps of the 
galaxy to date. 



2.2 CO observations 

The CO 1 — ► and 2 — > 1 observations of the inner ~ 8' x 2' 
region of the Circinus galaxy were performed with the SESTn at 
La Silla, Chile. In total, a grid of 7 x 27 points separated by 20" 
was observed along the major axis and in 7 x 13 points along the 
minor axis. The integration time was 15 minutes per grid point (in 
3 minute scans), giving a typical r.m.s. of T\ ^ 30 mK per 0.7 
MHz channel. On June 7-9, 1999, we observed the central three 
strips, at position offsets of -20", & +20" along both the major 
and minor axes and on April 22, 2000, we observed the remaining 
positions. The CO 1 — ► and 2^1 transitions were mapped 
simultaneously with the 115 & 230 GHz (IRAM) receivers, tuned 
to single-sideband mode. Typical system temperatures, on the T%- 
scale, were 300 - 400 K at 1 15 GHz and 300 K at 230 GHz. The 
back-ends were acousto-optical spectrometers with 1440 channels 
and a channel width of 1.8 krns - . We used dual-beam switching 
with a throw of about 12' in azimuth, and pointing errors, using the 
SiO maser W Hydra, were typically 3" r.m.s. on each axis. The half 
power beam-widths (HPBWs) are 45" for CO 1 -> and 22" for 
CO 2 — »■ 1. The intensity was calibrated using the chopper- wheel 
method. For all of the observing runs the weather was excellent, 
and only the removal of linear baselines, using the XS package, 
was required. 



1 The Australia Telescope is funded by the Commonwealth of Australia 
for operations as a National Facility managed by CSIRO. 

2 Operated until 2003 by the European Southern Observatory (ESO) and 
the Swedish National Facility for Radio Astronomy, Onsala Space Obser- 
vatory, Chalmers University of Technology. 
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Figure 1. High-resolution Hi moment maps of the Circinus galaxy. Left: The Hi distribution; the contour levels are 0.4, 1, 2, 3, 4, 6, 8, 10, 12, 18 and 20 
Jy kms _1 beam _1 . Right: The mean Hi velocity field; the contour levels range from 300 kms -1 (NE) to 570 kms -1 (SW) in steps of 10 kms -1 . The 
synthesised beam (124" X 107") is shown in the bottom left corner. 
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Figure 2. Low-resolution H I moment maps of the Circinus galaxy as obtained from the deep Parkes H I survey of the Zone of Avoidance IJuraszek et alj 
2000; Henning et al., in prep.). Left: The Hi distribution; the contour levels are (0.5, 1, 2, 5, 10, 20, 30, 40, 50, 60 and 66) X 7 Jy kms -1 beam -1 , where 
7 Jy kms -1 beam" 1 correspond to an Hi column density of 10 19 cm -2 . The contour levels were chosen to match those by Freeman et al. (1977; their 
figure. 6) above their detection limit of TVhi = 5 x 10 19 cm -2 ; our detection limit is a factor ~ 10 lower. Right: The mean H I velocity field (masked at TYjji 
= 5 X 10 18 cm -2 ); the contour levels range from 324 kms -1 (NE) to 534 kms -1 (SW) in steps of 10 kms -1 . The cross marks the center of the Circinus 
galaxy as given in Table[T] The gridded beam of 15 f 5 is shown in the bottom left corner. 
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3 RESULTS 

3.1 Gas distributions 

3.1.1 The H\ maps 

Fig. Q] shows the mosaiced, high-resolution H I distribution and 
mean H I velocity field of the Circinus galaxy as obtained with the 
ATCA. The equivalent low-resolution Hi images obtained from 
the HIZOA survey are shown for comparison (Fig. [2}- Our maps 
clearly show the enormous extent of Circinus, while exhibiting the 
asymmetric outer envelope, radial position an gle change and th e 
inner ~30' of the single-pointing H I image (see Jones et al J 19991) . 
2MASS JHK-band images and photometry Jjarrett et aJ1l2003h 
show an infrared extent of ~ 16' for Circinus, similar to the opti- 
cal Holmberg radius as extrapolated by Freeman et al. (1977) [see 
TableQ]. The H I diameter is at least a factor of five larger than the 

largest estimated optical and infr ared diameters of Circinus. 

As noted in previous maps iFreeman et al.ll 19771 ; Ijones et al.l 

1 19991) . the ellipticity and position angle changes with radius are 
indicative of a gentle warp. We also note some asymmetries in the 
gas distribution and velocity field with the extended H I emission 
towards the SWWfj being prominent. This pronounced emission in 
the more tenuous gas is in the same direction as the extension in 
the CO emission (discussed in Sect. l3.3.Tb . 



3.1.2 The CO maps 

In Figs. [3] and [4] we show CO contour maps of the fully mapped 
region, where we clearly detect CO emission out to +140" (2.9 
kpc) along the NE (I A = 4.5 ± 1.2 K kms -1 ) and -160" (3.3 
kpc) along the SW (I a = 3.8 ± 1.1 K kms -1 ). Note that, like 
lElmouttie et al.l ( 119980 , we detect very faint emission further out 
along the axes, but since these are detected at < 2a (e.g. Ia — 

2.1 ± 1.0 K kms -1 at +180" along the major axis), they are not 
considered further. 

Along the minor axis, CO is clearly detected to 80" along 
both the NW and SE (I A = 4.8 ± 1.3 K kms" 1 and I A = 

4.2 ± 1.2 K kms -1 , respectively). This deprojects to values of 
3.3-6.7 kpc for i = 65° to 78° (Table [j}. The lower end of this 
range agrees well with the value obtained from the major axis giv- 
ing i = 60° ± 6°. Note, however, that within the central 5 kpc, the 
H I velo city field displays an elongated structure and non-circular 
motions Uones et al J 19991) . This would have the effect of lowering 
our estimates of the inclination of the "disk", although the molec- 
ular ring+outflow model supports the high inclinations within the 
central kpc, with the velocity field within this region not exhibiting 
the same strong twists as on the larger scalqj. 

CO 2 — > 1 is detected at similar main-beam temperatures as 
the 1 — > transition at ~3.2 kpc and ftsl.l kpc along the major 
and minor axis, respectively, giving i = 65° ± 8°. These values 
suggest that the large-scale CO shares a similar inclination to the 
H I and o ptical disk (Table[T|>, and considering the value of i = 78° 
found by Curran et al. ( 1999), indicates that the CO distribution, 
like the Hi dJones et al . 1999), is warped. 



3 This is more apparent in the colour version of moment zero map, avail- 
able from http://www.phys.unsw.edu.au/~sjc/circinus/ 

4 Due to the low resolution, our first moment CO maps are fairly feature- 
less and the reader is referred to the over-sampled, higher sensitivity CO 
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Figure 3. The CO 1 — + map of Circinus. The contour levels are up to 80 
K kms -1 in 10% steps. The SEST antenna gain is 27 JyK -1 at 115 GHz. 
In this and Fig.|4]the position angle is 214° , and the dots mark the observed 
positions. 



3.2 Gas masses 

3.2.1 Atomic gas mass 

In Fig. [3] we show the global integrated H I spectrum of Circi- 
nus obtained from our mosaic. The integrated H I flux density of 
Fhi ~ 1300 Jy compare s with 1000 Jy kms -1 from the pri- 
mary beam limited image dJones et al.ll 19991) and Fm =1870 and 
1960 Jy kms" 1 from the single dish H IZSS and HIZOA maps 
JHenning et alJuOOOuJuraszek et alj20001) . respectively. Our value 
gives a mass of M H i = 5.5 ± 0.2 x 10 9 M R , using M m - 
2.36 x 10° D 2 F m M© dWildl 1 9521: iRobertsll 1 962J) and D = 4.2 
Mpc. 

Like the previous results, the spectrum shows a prominent 
double-horn shape, with more emission from the fast moving (at 
v « 330 kms -1 in Fig. |5J gas in the approaching (NE) segment. 
The emission peaks occur at 332 and 540 kms -1 , giving ±104 



1 map o flCurran et aT] Jl999l) . 



the sy stemic velocity of 439 ± 2 km s -1 derived bv lFreeman et al.l 
1 19771 [see also Sect. 13.3.31 and indicates that the asymmetry in the 
profile does not dominate the rotational dynamics of the gas. 
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Figure 4. The CO 2 — » 1 map of Circinus. The contour levels are up to 78 
K km s~ 1 in 10% steps. The SEST antenna gain is 40 Jy K~ 1 at 230 GHz. 
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Figure 5. Integrated H I spectrum of the Circinus galaxy. The velocity res- 
olution is 4 kms -1 and the integrated flux density is Fjji = 1330 ± 50 
Jy kms -1 . All velocities quoted throughout this paper are heliocentric. 
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Figure 6. Integrated CO 1 — > (black) and 2 — ► 1 (coloured) spectra of 
Circinus. The velocity resolution is 10 kms -1 and the antenna tempera- 
tures have been corrected by main beam efficiencies of 7? m b = 0.7 and 
0.5, respectively. The Galactic line close to km s — 1 occurs in th e western 
part of the map clo se to known Ga lactic emission IDame et al.lll987l) . as 
previously noted bv lElmouttie et al.1 1 1998). 



3.2.2 Molecular gas mass 

Fig.[6]shows the integrated CO 1 — ► and 2 — * 1 spectra summed 
over the full map. For both transitions we see the excess of low 
velocity gas in the receding peak, as was evident in the molec- 
ular ring model of ICurranet alj J 19981) . This contravenes the Hi 
spectrum (Fig. [5}. indicating a significantly larger molecular frac- 
tion in the SW portion of the ring. For CO 1 — > the average 
velocity integrated line intensity is I a — 6.28 ± 0.07 K kms - . 
Converting this to main beam intensity, via ?7 m b = 0.7 ± 10% 
at 115 GHz, gives an average of J m b = 9 ± 1 K kms -1 over 
the whole map (Fig. resulting in a CO 1 — » luminosity of 
Leo 1-.0 = 400 ± 40 K km s _1 kpc 2 , where we have added an ad- 
ditional half beam-width (S m b) to each axis to ensure we account 
for all of the mapped regiorQ 

Applying the Galactic N iio/Ico conversio n ratio of 2.3 ± 
0.3 x 10 20 cm~ 2 /[K kms" 1 ! Jstrong et al.ll988i) and a 36% mass 
correction for helium (so that 1 M = 4.38 x 10 56 H2 molecules), 
gives a total molecular mass of Mh, = 2.0 ± 0.5 x 10 9 M Q . This 
is the same as the value found by Elmouttie et al.1 jl99oti 6 l and sim- 
ilar to that obtained bv lCurran et alj dl998l) over the central region 
only. This therefore suggests that the previous ±60" x ±60" maps 
contained most of the CO emission. 

Since this mass is comparable with t he dynamical mass , 
M dyn = 3.3 ± 0.3 x 10 9 M within 560 pc. lCurran etalU 1998b 
raised the possibility that the Galactic conversion ratio was in- 
applicable in this case, as had been previously noted for other 
galaxies (e.g. iMalonev & B lack 1988; Malonev 1990; Shier et al. 



5 For small angles the fractional uncertainty in this conversion is approxi- 
mately equal to that of the distance estimate, which is =3. = £4^ = 0.19 
(Sect. [T}. This results in fractional errors of =rf- ~ 0.4 for the gas mass 
estimates and m 0.2 for the dynamical mass estimates, although these un- 
certainties will cancel somewhat when comparing these masses (i.e. if one 
is overestimated so is the other). 

6 They actually quote Mh 2 = 1.1 ± 0.2 x 10 9 M , giving 
MH 2 /Mdyn = 0.05 ± 0.01, using a -/Vh 2 //co conversion ratio of 
1 x 10 20 cm- 2 /[K kms" 1 ]. 
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Figure 7. Overlay of the global (black), ring region (red) and central 
(dashed blue) CO 1 — > spectra. All have been corrected by 27 Jy K — 1 
and the ring region and central spectra scaled up by factors of 1.5 and 10, 
respectively. 



Il994l ; |Pahmer]|i995l ; iGenzel et alj|l995h . In order to check that 
this apparently high molecular mass was not the result of integrat- 
ing the luminosity over r > 560 pc, we averaged the intensities 
over the central 25 positions of Fig. [3] i.e. to ±40" x ±40", 
which gave J m b = 48 ± 5 K kms -1 . Deconvolving O m b from 
this ±40" gives a radius of 680 pc in which a molecular mass of 
Mh 2 = 1.0 ± 0.2 x 10 9 Mq is found, thus either suggesting a gas 
mass fraction of « 30%, within this radius, or a conversion ratio of 
Nn 2 /Ico < 1 x 10 20 cm~ 2 /[K km s" 1 ]. 

Note that lElmouttie et al.l d 19981) previously obtained a much 
lower molecular gas mass fraction from the central CO 1 — » 
spectrurrM Estimating the dynamical mass over the central beam 
in a similar fashion for our data, i.e. assuming that the maxi- 
mum rotational velocity occurs at the edge of the beanQ gives 
Mdyn ~ 3.1 x 10 9 M Q within the central < 440 pc, compared to 
the « 3.9 x 10 9 M Q found bv lElmouttie et al] Jl998h - The molecu- 
lar mass over th e same region is Mh 2 ~ 4.4 x 10 s Mq, giving, like 
lElmouttie et al.l d 19981) , Mn 2 /M dyn as 0.1. The large increase in 
the gas mass beyond this radius (previous paragraph), increases the 
gas mass fraction by a factor of three (to A/h 2 /M dyn fa 0.3) within 
680 pc, thus being consistent with the presence of a molecular ring 
in the central regions of Circinus. Note, however, that an overesti- 
mate in the dynamical mass is possible when there are non-circular 
motions present, due to the presence of non-a xisymmetric potentia l 
(e.g. a bar), particulalry in the central regions (Koda & Wada 2002). 
Therefore, the ring interpretation should only be considered as one 
possibility. Fin ally, we see that none of the central 1 — > 0, 2 — > 1 or 
3^2 spectra dCurran et al. 1998, 2001) exhibit the horns apparent 
in the global spectra. This difference is also seen in our data, where 
in Fig [7] we overlay the global CO 1 — ► spectrum upon those of 
the centre position (within the central 45" beam) and the molecular 
ring region (within the central ~ ±40 x 40" mapped positions). 
This shows that our large-scale CO maps detect the more extended 
(;£ 1 kpc) gas, which rotates at ~ ±130 kms -1 (see Sects. [3331 
andl4~il. 

Converting the average CO 2 — »• 1 velocity integrated line in- 



tensity of /a = 3.82±0.04Kkms" 1 ,via77 mb = 0.5±20%at230 
GHz, gives 7 m b = 7.7 ± 1.5 K kms -1 resulting in a CO 2 — ► 1 
luminosity of Leo 2->i = 350 ± 70 K k ms -1 kpc 2 . From pre- 
vious observations of the central position jjohansson et al.ll 1991c 
ICurran et alJuOOll) . ^o iHo intensity ratios of close to unity were 
noted and Figs.[6]and|7]indicate that this ratio drops slightly in the 
horns of the CO profile, which traces the extended gas beyond the 
ring. 

The H 2 (derived from CO 1 -> & 2 -» 1) and Hi column 
densities along the major axis are shown in Fig.[8] Since the former 
is derived from both the CO 1 — > and 2 — > 1 transitions using 
the same 7Vh 2 /Ico conversion ratio, it represents a profile of the 

CO 2->l 
CO l->0 



intensity ratio along the major axis. Over the range where 



we have clear detections (for H2 column densities of 7Vh 2 ^ 10 21 
cm -2 ), we see that the high ^q ^q ratio holds. This may suggest 
that the conditions (T kin ~ 50 - T50 K, N CO ~ 0.2 - 20 x 10 18 
and r t H, ■> 10 3 cm" 3 ) found ove r the central beam dCurran et al.l 
1998, 2001 : lffitschfeld et al.ll2008h could hold over the bulk of the 
molecular gas, although these are high temperatures to maintain 
over a 3 kpc radius. Within this region we see very high molecular 
fractions of / = 27Vh 2 /[27Vh 2 ± N m ] ~ 1, the value of which 
drops to / < 0.7 beyond the limit of the clear CO detection range 
(> 3 kpc). That is, in the case of Circinus we witness the limiting 
atomic gas column density of Nm S, 10 21 cm~ 2 , which is typi- 



7 Figure 9 of ICurran et alj 09981) shows that this is a fair assumption in 
this particular case. 



cal fo r disk galaxies (e.g. ICavatte et alJI 19941 ; iRhee & van Albadal 
1996), with gas above these column densities being converted into 
the m olecular state under the favourable cond itions close to the nu- 
cleus JYoung & Scovillelll99ll;ISchavell200il) . 



3.3 Gas dynamics 

3,3.1 Position angle variations 

Previously, Ijones et alj dl99 9) noted two large spiral arms origi- 
nating at the end of the atomic ga s bar at a radius of ~ 5 kpc 
and, as seen in the H I distribution d-Iones et al.ll 19991 and Fig. [TJ, 
the bar itself is enclosed by a large ring of atomic gas. Referring 
to the CO emission (Figs. [3] & |4J, being limited to < ±3' the 
molecular gas extends to only the very central regions of the H I 
(Fig. [U left). However, similar to the Hi (Sect. l3~l~Tt . there ap- 
pears t o be an extension in th e more tenuous emission towards the 
SWW. Elm outtie et alj d 19981) also noted this "disturbance" in the 
kpc-scale CO emission, but since this was due to a single position 
spectrum, they deemed it doubtful. We find, however, that feature is 
visible over several spectra and at the peak has an integrated inten- 
sity of 7a = 20.3 ± 1.3 K kms - . This occurs at a distance of 1.8 
kpc from the centre position with a position angle of PA ~ 236°, 
thus exhibiting a "warp" from the ~ 214° position angle of the 
molecular ring (Fig. [9}. Su ch a change in pos ition angle is also seen 
in the atomic gas (see also ljones et al.l ll999 and Sect, l3,33t on the 
10-kpc scale and Fig. [9] illustrates that similar changes also occur 
at the inner edge of the accretion disk. Although contained within 
the central 0. 1 pc, the accretion disk shares a similar position an- 
gle (209 ± 3°) to the molecular and atomic rings, while the disk 
at 0.4 pc (236 ±6°) has the same position angle as the extended 
molecular (> 1 kpc) and atomic gases (> 10 kpc). 

These abrupt changes in position angle could be due to 
large non-circular motions, perhaps caused by bar potentials in 
the inner parts of the galaxy w here the CO resides , as has been 
found for many other spir als dScoville et alj I 1993k ]Regan et al. 
1999; Sakamoto et~al]|l999l) . including the Milky Way dFuxHl999l 
iDame et alj 1200 II) . In these cases, the molecular gas is found to 
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Figure 8. H I (crosses) and H2 (circles) column densities as a function of galactocentric radius (NE to SW along the major axis). The H I column densities 
are calculated in the optically thin regime from a brightness temperature of 45.93 K Jy -1 per the main-beam solid angle of 3.52 X 10 — 7 sr. The H2 column 
densities have been calculated using Nn 2 /Ico = 2.3 X 10 20 cm - 2 /[K km s ~ 1 1 (Strong et al . 1988): For CO 2 — + 1 (unfilled circles), we use a main-beam 
efficiency of rj m b = 0.5 and for CO 1 — » (filled circles), r? m b = 0.7. The dashed lines represent the offsets of the last clear (> 3<r) CO detections. 
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Figure 9. The "zig-zagging" of the positio n angle on various sca les in Circinus. Left: The molecular gas - the inset shows the H2O maser emission from the 
warped accretion disk in the inner parsec (Greenhill et al. 2003), which is superimposed upon the deconvolved (to a resolution of 10" using the method of 
Rydbeck 2008) map of the inner CO 2 — ► 1 emission (Curran 1998). This resolves the ring in the NE-SW diagonal, as well as the perpendicular outflow 
iCurran et al . 1998, 1999). Middle: The CO 2 — > 1 first order moment map (due to binning of the data, the weaker extreme velocity channels are lost, see Sect. 
13.3.3) . Right: The central H I first order moment map. In the latter two maps we show the beam in the bottom left corner, 22" and 124" X 107", respectively. 
In the colour version of the figure, blue is approaching, red is receding and green is systemic gas. 



abruptly change direction when transversing the dust lanes o n the 
way to the nuclear ring dRegan et alj |l999; Dame et al feOOll). and 



in the case of Circinus, a dusty nuclear bar JMarconi et al 



199. 



nestling within the molecular ring (see figure 1.7 of Curran 200(T 



has been found. For Circinus it is apparent that a ring forms at 
the ends of each bar and a tendency for strong bars to terminate 
in rings has been noted in the dust structure of 75 galaxies by 



IPeeples &M artini (2006). Furthermore, as is seen for other galax- 
ies (e.g. Sakamoto et al. 1999) and from numerical simulations 
of barred potentials (Wada 1994; Bureau & Athanassoula 1999; 
lAthanassoula&Bureaulll999l ; lKoda&Sofuell2006l) , the nuclear 
bar may provide a means for transporting the gas to the nucleus 
of Circinus. Further support is lent by the gas mass fraction of 
MH 2 /A^dyn ~ 0.1 within ~ 400 pc, derived above, which is con- 



sistent _wMiJhe_bar_driven transport of molecular gas to the nu- 
cleus (Sak amoto et alj|l999l) . For Circinus, we therefore suggest, 
that upon each collision with a bar or ring, the gas undergo a jolt 
which shifts the position angle by ~ 20°. As illustrated in Fig. [9] 
these same twists and turns in position angle occur on wildly dif- 
ferent scales and this may be the first case that the orientations and 
kinks are seen to be repeated all of the way down to the sub-pc 
accretion disk. 



http://nedwww.ipac.caltech.edu/level5/Curran/frames.html 
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3.3.2 Position-velocity diagrams 

In Fig. [10] we show the Hi and CO 1 — * position-velocity (p-v) 
diagrams of Circinus along its major axip As e xpected, the H I 
distribution is similar to that of I Jones et al.l J199 9). with a rising 
component extending to a radius of ~ 400" (s a 8 kpc), be i ng fol - 
lowed by a flatter component. We also see, like ljones et alj (1999). 
that along the major axis the gas is more extended in the NE than 
the SW (also visible in Fig.[8j, being clearly visible to 2000" (« 40 
kpc) and 1600"(~ 30 kpc) respectively. 

From the CO distribution it is apparent that the molecular gas 
traces a different kinematical feature to the bulk of the atomic gas, 
in which there is a significantly higher velocity gradient and a larger 
velocity range, giving the difference in the global H I and CO pro- 
file widths (Fig. II It. High velocity wings in the CO distribution of 
other galaxies have previously been attributed to high dispersion or 
rapid rotation of the molecular gas in the inner disk dSofudll992l ; 
iTutui & Sofudll999l) . Zooming in on the H I emission profile over 
the CO emission region (Fig. I lit, we see that this accounts for the 
weak large velocity wings of the global H I emission with a width 
close to that of the CO profiles. We therefore interpret this emis- 
sion and the wings of the CO profiles to t he fastest part of th e ring 
(> 130 km s" 1 ), within the inner 600 pc dCurran et all 19981) . 

The CO p-v d istribution is very similar to that of 
lElmouttie et al.l 1 119981) : Along the major axis we also observe 
'rigid' rotation within the central ±60", as well as the flat rotation 
at distances beyond this. The rigid rotation is consistent with the 
model of the sub-kpc molecular ring, where the velocity increases 
to a d eprojected maximum of 180 kms -1 at 400 pc dCurran et al.l 
1998). When convolved with the 44" beam, this transforms to an 
observed offset of «50", thus suggesting that the rigidly rotating 
inner ±60" corresponds to the molecular ring. Finally, the loca- 
tions of the peak emission at velocities of 294 ± 2 kms -1 and 
575 ± 5 kms -1 , give symmetrical rotation about 435 ± 5 kms~ , 
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- dashed orange) emission profiles. The spectrum of the H I over the CO 
emission region has been scaled up by a factor of 25. 



in excellent agreement with the systemic velocity derived from the 
Hi (Sects. l3Xn and l333l 

Beyond the inner ±60", the flat rotation along the major axis 
is seen to occur at ?»300 kins" 1 (NE) and r;560 kms -1 (SW). 
Comparing these to the systemic velocity, gives a projected ve- 
locity of about ±130 kms -1 , a value similar to that found by 
lElmouttie et al.l (JJ998). This corresponds to a deprojected veloc- 
ity of «130 to 150 kms -1 (for i = 78° to 60°. Sect. l3~L2l , over a 
deconvolved distance of ~800 pc to 4 kpc. The inner edge of this 
large-scale, solid body molecular disk matches well with the inner 
rigidly rotating ring which has velocities of < 150 kms -1 at >600 
pc dCurran et alll99cj> . 



9 Due to the lower resolution and sensitivity to the inner regions, the dia- 
gram for the m inor axis is fairly fe atureless and the reader is referred to the 
p-v 



diagram o flJones et alj (1999). 



3.3.3 Til ted- ring model 

In order to model the velocity field in Circinus a tilted ring model 
was fitted to the first order moment maps. This was performed 
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as per the modelli ng of iRogstad et al.l \l9l4) : iBegemarj Jl989l) ; 
iThean et al.l dl997l) . using with the MPS task gal, where by as- 
suming that the velocity field was constant over each 120" an- 
nulus (~ lx the beam-size), we fitted concentric annuli. In the 
first iteration, all six of the fitting parameters were allowed to vary; 
the central position coordinates, the position angle of the reced- 
ing major axis, the inclination of the orbital plane, the systemic 
velocity and the maximum rotation velocity. The fits to the inner 
four annuli, from 120" to 600", where the velocity field is most 
symmetric, each gave the same central position and systemic ve- 
locity of ajaooo = 14 h 13 m 10 s , 5j2ooo = -65° 20' 49" and 
whi = 440 kms -1 , respectively. That is, within 30" (a quarter 
of the HP BW) and 1.5 km s" (half a channel) of the values ob- 
tained by Ijones et alj q 19990 . Adopting the Jones et al. values to 
drive the tilted ring model, we then ran another iteration of gal, 
keeping these values fixed and looping round the annuli in order to 
fit the full velocity field. This was done out to 2460", as beyond 
this the data proved too poor to fit a model^j. This gave a three 
parameter fit very similar to that of the initial model for the whole 
H I velocity field, as well as separate fits for the approaching and 
receding sides. The results are shown in Fig. 1 1 3l L=L L 

To model the molecular gas, datacubes of the CO 1 — ► and 
2 — » 1 were constructed in CLASS, in first order moment maps with 
5-point hanning along the spectral axis. These were smoothed in 
order to improve signal-to-noise ratio to a point where the data were 
useful. We then performed the modelling as for the Hi, fitting a 
rotation curve with the constant velocity model and using the centre 
position and systemic velocity derived from the Hi data. Again, 
we used annuli of lxHPBW width (45" for CO 1 -> & 22 " 
for 2 — > 1), out to 200" and 170", respectively, where the data 
appeared reasonable, although none of the CO first moment maps 
where of nearly as high a quality as for the H I. 

This is the first time that a tilted ring model has been applied to 
the molecular gas in this galaxy and we see that for the inner rota- 
tion curve (Fig. [12] top), the fit is consistent with the Hi data, with 
the outer values of ~ 130 km s -1 agreeing well with the those de- 
rived by assuming circular orbits inclined at 2; 70° (Sect. [3~3~2t . 
Being restricted to radii of 2 900 pc, however, the model can- 
not trace the high velocities of the inner molecular ring discussed 
above. 

The CO model fits to the position angle are more awry, al- 
though still within < 20° of the H I model (Fig. [12] middle), but 
larger than the changes in position angle discussed in Sect. 13.3.11 
For the inclination (Fig. [12] bottom), some CO fits are reasonably 
consistent with the Hi model within the inner 160", but jump up 
to ~ 100° beyond this, at radii where the drop in position angle is 
seen. We therefore attribute this to the poorer first moment maps of 
the CO emission, which were more difficult than the H I to model, 
although, again, we do expect the inclination to rise within the inner 
kpc (Table[2j. 

Comparing the H I model with that of Ijones et alj dl999l) , we 
see that out to a radius of 1700", where their model was deemed re- 
liable, there is excellent agreement in all three parameters, with the 
rotation curve, as well as all of the peaks and troughs in the position 
and inclination angles, being reproduced. Interestingly, the offsets 
of the receding model fit from those of the approaching half and the 



10 Due to the blanking of some "hot" pixels in the receding side this is only 
modelled to 1800". 

11 The model velocity field and the residuals are available from 
http://www.phys.unsw.edu.au/~sjc/circinus/ 
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Figure 12. The three parameter fit to the H I (filled circles) and CO (trian- 
gles) velocity fields. The small circles show the model fit to the approach- 
ing side of the H I, the crosses the receding side, the filled triangles the CO 
1 -» fit and the unfilled triangles the CO 2 — ► 1 fit. 



Table 2. Inclination angles, i, at various radii in Circinus. 



Species 



Radius 



Reference 



H 2 

CO 



Ho 
Hi 



0. 1-0.4 pc 

0.1-0.6 kpc 

3 kpc 

3 kpc 

4 kpc 
8 kpc 

26 kpc 
39 kpc 
39 kpc 
47 kpc 



edge-on 

78° ± 1° 
> 63° ± 5° 

60° ± 6° 

40° ± 10° 
ft 70° 
ft 40° 
ft: 60° 

65° ± 2° 
^ 50° 



Greenhill et al. (2003) 
Curran et al. (1998) 
Elmouttie et al. (1998) a 
This paper b 
Elmouttie etal. (1998) 
Jones et al. (1999) 
Jones et al. (1999) 
Jones et al. (1999) 
Freeman et al. (1977) 
This paper 



Notes: "Elmouttie et alj 11998) actually quote an inclination angle of at 
least 73°, however we obtain 63° based on the CO extent along the ma- 
jor and minor axes quoted. By compari ng the scale size s of the minor and 



major axes emission and applying, like lElmouttie et alj 419981) . an uncer- 
tainty of one mapped position per each full axis (see Sect. 13.1.21 ) 
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Figure 13. The rotation curve over the major axis of the Circinus galaxy. 
The dotted curve shows the Keplerian fit to the H2O masers I Greenhill et al. 
2003), the stars and squares sho w the stars and Nil emission line data 
iMaiolino et al.|[l99 8; Oliva et al] 1 19941) , respectively, and the circles the 
show the H2 data of Da vies et alj 419981) . The filled blue symbols represent 
the NE (approaching) data and the unfilled red symbols the SW (receding) 
data over these regions. Farther out, th e symbols are as pe r Fig. 1 121 the black 
triangles represent the CO ring model ICurran et alj 19981) . the coloured tri- 
angles the larger-scale CO (1 — ► — filled green & 2 — > 1 - unfilled 
orange) and the small circles the H I. 



whole H I velocity field are also reproduced over ~ 1100 — 1800", 
confirming the as ymmetry of the outer Circinus disk as found by 
Jones et al. ( 1999). Beyond the inner 1700" (which is discussed in 
detail bv ljones et al. 1999), both the position and inclination angles 
settle close to the outer values previously found. With regards to the 
rotation curve, after one more small climb, this remains fairly flat, 
indicating that the fastest rotating atomic gas is associated with the 
large-scale diffuse emission and not within the sub-kpc ring (Sect. 

I335t , 



4 MASSES IN CIRCINUS 

4.1 Dynamical masses 

In Fig.Q~3]we combine the rotation curves from the tilted ring mod- 
els with those from the literature to show the dynamics over the 
whole galaxy. The overall profile is fairly typical for spirals (type 
Sa - Sc) with a steep initial rise within the inner few kpc, beyond 
which the velocity plateaus JRubin et alJl98CLIl982Lll985tlBosm3 
ll981al lbT). 

Specifically, in the inner 0.4 pc we see Keplerian rotation in- 
dicative of dynamics which are d ominated by a compact massive 
object (see iGreenhill et alj |2003). Beyond this, from fa 10 pc to 
ss 400 kpc the velocity scales directly with distanca 12 L despite 



12 There is some scatter, particularly at the lower radii, due to the ef- 
fects of seeing and in Figs. |13| and ll4l we have therefore excluded the cen- 
tral three points of iDavies et alj |l998l) Iwhere r< 1", cf. a seeing of 
>1.4"] and the central point of iMaiolino et alj (1998) [where the seeing 
is f»0.5" - 0.6"]. As well as the seeing, uncertainties in the velocities are 
quoted: 9 to 15 kms for the stars and 30 kins" 1 , for the Nil emission 
IMaiolino et alj 19981) , which we do not show for the sake of clarity. While 
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Figure 14. The dynamical mass distribution of the Circinus galaxy. The 
symbols are as per Fig.[l3] This confirms the underestimate in the velocities 
from the CO 2 — > 1 tilted ring model, although the outer three 1 — ► 
estimates appear quite reasonable. 



both the H2 and C O intensity dropping si gnificantly over the in- 
ner few hundred pc I IDavies et al.ll 19981 and lCurran et alj|l998l . re- 
spectively), probably due to consumption by the starburst activit y 
within the central kpc JMarconi et alj [l994; Elm outtie et al.ll 19981) . 
The outer edge of the starburst/molecular ring is marked by 
a sharp drop in the rotational speeds, although probably not to the 
extent as indicated by the CO 2 — ► 1 tilted ring modeo, which is 
confirmed by the dynamical mass distribution (Fig.ll4t. Beyond the 
molecular ring, velocities settle at ~ 150 kms -1 , where the gas 
follows the general, large-scale solid body rotation typical of disk 
galaxies, although Fig. [14] may suggest that there is not much ad- 
ditional (gravitating) matter between the ring and the large-scale 
atomic gas. Beyond this plateau, the mass of the galaxy clearly 
keeps climbing beyond the outer edge of the disk, although with 
less of a gradient than the inner kpc, confirming the differences 
in the p-v distributions between the atomic and molecular gases 
(Fig. UOt. Visually, there is a drop of 3 magnitudes within the cen- 
tral 1 kpc and with a Holmberg radius of 10 kpc (Freeman et al. 
1 19771) [see also Sect. |3.1.T| , it is clear that the halo must be domi- 
nating the outer dynamics of Circinus. 



4.2 Gas mass fractions 

In Table [3] we summarise the various gas a nd dynamical masses 
found for the Circinus galaxy. ICurran et alj i 19981) reported a gas 
mass fraction of Mn 2 /M<iy n ~ 0.5 within the central 560 pc, 
which called into question the Galactic J Vh, /Ico conversio n ra- 
tio of « 2.3 x 10 20 cm~ 2 /[K kms -1 ] dstrong et al.lll988h . We 
see, however, that this does not apply ov er the whole of the galac- 
tic disk, falling to the nominal < 0.1 JMihalas & Binnevlll98ll) 
beyond the molecular ring, and upon reaching the Holmberg ra- 
dius, the gas mass fraction fall s to « 0.03. Combining the atomic 
gas mass Jjuraszek et alj|20001) with the molecular gas mass gives 



these data are characterised by linear fits on linear plots, in the log plot 
uncertainties in the position are hugely exaggerated. 

13 Also, the decrease would be much more gentle with radius than evident 
from the logarithmic scale of the plot. 
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Table 3. Gas and dynamical masses at various radii in Circinus. 



Mass [M ] 


Radius 


Type 


Reference 


Ri 4 X 10 8 
< 4.7 x 10 9 


< 440 pc 

< 440 pc 


H 2 a 
Total b 


Elmouttie et al. (1998) 
Elmouttie et al. ( 1998) 


1.0 ±0.2 x 10 9 
3.3 ± 0.3 x 10 9 


< 680 pc 

< 560 pc 


H 2 a 
Total 


This paper 

Curran et al. (1998) 


2.0 ±0.5 x 10 9 
si 2.1 x 10 10 


< 5 kpc 

< 5 kpc 


H 2 a 
Total 


This paper 
This paper 


1.9 ±0.1 x 10 9 
«6.2 x 10 10 


< 10 kpc 

< 12 kpc 


Hi 

Total 


This paper 
This paper 


^ 3.8 x 10 9 
1.9 ±0.3 x 10 11 


< 35 kpc 

< 35 kpc 


Hi 

Total 


Jones et al. (1999) 
Jones et al. (1999) 






7.2 ±0.5 x 10 9 

1.3 ±0.2 x 10 11 


< 40 kpc 

< 40 kpc 


Hi 

Total 


Freeman et al. (1977) 
Freeman et al. (1977) 






«8.1 x 10 9 
2.6 ±0.6 x 10 11 


< 60 kpc 

< 50 kpc 


Hi 
Total c 


Juraszek et al. (2000) 
This paper 



Notes: a Applying the N^ n / Iqo conversion ratio of 2.3 x 10 20 cm 2 /[K 
km s ~ * 1 fetrong et al J l 1 9881) Isee also Sect. 1331. b CO using i = 63° (Ta- 
ble[2) lElmouttie et alH1998t) quote < 3.9 ± 0.2 X 10 9 M for i = 73°]. 
c assuming that the uncertainty in the tilted ring model, which is not given 
by gal, is much less than the uncertainty due to the distance estimate of 
Circinus (see footnote|5j. 



Table 4. Mass inventory of the Circinus galaxy. The masses quoted are the 
(most) global values. 





Mass [ M Q ] 


Radius 


Molecular gas 
Atomic gas 
Ionised gas 
Stars 
Dust 


2.0 ±0.5 x 10 9 

rs 8.1 x 10 9 

~ 10 6 

-10 11 

~ 10 6 


< 5 kpc 

< 60 kpc 
<"o.3 kpc 

< 10 kpc 

< 60 kpc 


Total 


-10 11 


< 60 kpc 


Dynamical mass 


« 3 x 10 11 


< 50 kpc 



M g as = Mm + Mh 2 « 1 x 10 10 M Q , and so M gas /M dyn « 0.03 
holds over the inner ~ 50 kpc of the galaxy (Table |3), This insig- 
nificantly smaller than the value obtained by Elm outtie et al.l ( 1 998) 
[Mgas/Mdyn ~ 0.6], and is cons istent with the SAS3 (Hu bble type 
Sb, where M raB / M vn ~ 0.08.lYoung & Scoville l ll99ll) classifi- 
cation of Circinus dde Vaucouleurs et aljl99lh . 



4.3 Mass inventory 

In Table|4]we give the global values for the various masses in Circi- 
nus. The steep decline in column density beyond the central 3 kpc 
(Fig.[8jl, suggests that we have mapped most of the molecular gas, 
although from the gas-to-dust ratio, ISchmidt & Boiled j 19931) esti- 
mate Mh 2 « 0.6 x A/hi, which would giv e M H , » 5 xlO 9 M 
using the total atomic gas mass of I Juraszek et alj (2000). The 
ionised gas mass is obtained from the integrated Ha flux of the 
ionisa tion cone within the central 30 x 30 arcsec 2 JElmouttie et al.l 
1998). Since most of the ionised gas appears to reside in the ion- 
isation cone and star-burst ring, within a radius of « 300 pc 



JMarconi et al J 19941) . we do not envision the to tal ionised gas m ass 
competing with that of the neutral gas (recentlv. lRov et al J2008I es- 
timate an ionised gas mass of between 3 x 10 3 M© and 1 x 10 6 
Mq within a radius of 0.7 kpc). For the stellar mass content we have 
simply used the value typical for a large spiral and the dust mass 
is esti mated from the gas-to-dust ratio given by I Schmidt & Boiled 
(1993). From this, it is apparent that the stellar population could 
easily account for the observed dynamical mass, although, as dis- 
cussed in Sect. 14.11 the blue and visual brightnesses have all but 
petered out by time a radius of 10 kpc is reached. At this radius the 
dynamical mass is Mdyn ~ 5 x 10 10 M Q , with « 10 11 M Q being 
reached within the central 20 kpc CFig.l 14-b. 



5 SUMMARY 

We have studied the atomic and molecular emission in the Circinus 
galaxy through large-scale maps of the H I and CO 1 — > & 2 — > 1 
transitions. Our Hi mosaic represents the largest high resolution 
image of this galaxy to date, in which we recover 70% of the total 
H I flux, while resolving the detailed gas distribution. We have mo- 
saiced a field of ~ ±1° x ±1°, thus re covering m 3 0% more flux 
than the primary beam limited image of Uones et al.l J1999]) , giving 
the atomic gas distribution and dynamics to a radius of 50 kpc. The 
CO maps also represent the largest to date, spanning ±260" and 
±100" along the major and minor axes, respectively, at a sensitivity 
of « 10 mK per each 10 km s _1 channel. From our observations: 

(i) We confirm that the atomic gas in Circinus is distributed in 
a > 50 kpc irregular disk, which is a factor of > 5 times as ex- 
tended as the optical and infrared emission. A tilted ring model of 
the H I velocity field agrees very well with the know n warping in 
position angle and inclination, found by I Jones et al.l d 19991) , out to 
~ 30 kpc. Beyond this (to 50 kpc), the warping appears to dampen 
significantly. 

(ii) We obtain clear detections of CO at a galactocentric radius 
of 3 kpc. Although the position angles and inclinations yielded by 
the tilted ring model of the CO velocity field may be unreliable, 
by assuming circular orbits for the molecular gas we find that the 
molecular disk also undergoes a (« 20°) warp within the inner 
3 kpc, a radius at which it shares a similar inclination to the atomic 
gas. Note, however, that the apparent warp may be the result of non- 
circular motions, which are a disti nct possibility due to the presence 
of the atomic bar on these scales jjones et al.lll999h . 

(iii) The inner part of each gas structure also appears to maintain 
the same position angle on all scales; sub-pc (H 2 0), sub-kpc (CO) 
and < 10 kpc (H I). This may suggest a continuous transportation of 
material to the active galactic nucleus, while "kinks" are introduced 
when the outer diffuse gas encounters a bar. 

(iv) We confirm the finding of the ring model of ICurran et al.l 
(1998), that there is more molecular (and less atomic) gas in the 
receding (SW) half of the ring. 

(v) We find that the CO profile is wider than that of the H I in 
both the full-width half maximum and the total velocity spread. A 
comparison of the H I and CO position- velocity diagrams confirms 
that the molecular gas has different kinematics to the atomic gas. 
We believe that this indicates that the fastest rotating gas, which 
occurs within the sub-kpc ring, is indeed primarily molecular. 

(vi) We also investigate the global integrated atomic and molec- 
ular gas intensities, and find: 



(a) Circinus has a total molecular gas mass of Mh 2 



2 x 



12 S. J. Curran, B. S. Koribalski & I. Bains 







10 M<a, cf. the ato mic gas mass of Mm ~ 8 x 10 M 
jjuraszek et alfeood). 

(b) Previously, ICurran et al.l J1998I) found a gas mass frac- 



tion of Af gas « Mn 2 ~ 0.5Mdyn, within the central kpc, 



which cast into doubt the validity of applying the standard 
Galactic conversion ratio o f Nh 2 /Ico ~ 2 x 10 20 cm _2 /[K 
kms -1 ] dStrong et al.lll988D . However, over the global molecu- 
lar gas emission we obtain a more canonical (A/ gas < 0.1Md yn , 
iMihalas & Binnevlll981l) ratio of fa 0.03. This is much smaller 
than the value (M gas /Afd yll ~ 0.6) obtained bv lElmouttie et al.l 
J1998h . 

(c) Within the inner ~3 kpc, the gas remains predominantly 
molecular, with the molecular hydrogen column density ranging 



-,22 

10 21 



hydrogen column density of Nm < 10 21 cm -2 . Furthermore, 
the CO 2 -> 1/1 -> intensity ratio of ~1 previously observed 
over the central beam is also seen over the whole of this region. 

(vii) From the rotational dynamics of the gas, the steepest gra- 
dient of enclosed mass, dM/dr, occurs from wlO - 400 pc, al- 
though there is a low abundance of gas in this region (again, 
Afgas < O.lMdyn), probably due to consumption by the starburst. 
At > 400 pc, which is approximately coincident with the molecu- 
lar ring (Mgas/A/dyn ~ 0.5), there is a steep drop in v elocity. This 
continues to ~ 1 kpc, the extension of the atomic ring jjones et al.l 
1999). Beyond these radii, velocities climb once more, until set- 
tling into solid body rotation at ~ 150 kms -1 . This continues to 
the 50 kpc limit of our tilted ring model, although most of the blue, 
visible and infrared radiation emanates from the central 10 kpc, 
suggesting that the dynamics beyond this region are halo domi- 
nated. 
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